Positronium atom scattering by H2 in a 
coupled-channel framework 

P K Biswas^ and Sadhan K Adhikari + 
$ Departamento de Fisica-IGCE, Universidade Estadual Paulista 
13500-970 Rio Claro, Sao Paulo, Brasil 
+ Instituto de Fisica Teorica, Universidade Estadual Paulista, 
01405-900 Sao Paulo, Sao Paulo, Brasil 

February 2, 2008 

Abstract 

The scattering of ortho positronium (Ps) by H2 has been investigated using a three- 
Ps-state [Ps(ls,2s,2p)H2(X coupled-channel model and using Born approximation 

for higher excitations and ionization of Ps and B X S+ and b 3 S+ excitations of H2. We 
employ a recently proposed time-reversal-symmetric nonlocal electron-exchange model 
potential. We present a calculational scheme for solving the body-frame fixed-nuclei 
coupled-channel scattering equations for PS-H2, which simplifies the numerical solution 
technique considerably. Ps ionization is found to have the leading contribution to target- 
elastic and all target-inelastic processes. The total cross sections at low and medium 
energies are in good agreement with experiment. 
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1 Introduction 



With the availability of improved monoenergetic ortho positronium (Ps) beam, low-energy 
collision of exotic ortho Ps atom with neutral atoms and molecules is of interest in both physics 
and chemistry due to its vast applicational potential Recently, measurements of total 
ortho-Ps scattering cross section from various atomic and molecular targets (He, Ne, Ar, H 2 , 
N 2 , C4H10, and C5H12) have been carried out at low and medium to high energies @, ||, [|. 
Experimental investigations are also in progress for the measurement of pick-off annihilation of 
ortho Ps from closed shell atoms |J. Because of the composite nature and high polarizability 
of Ps, a reliable theoretical description of Ps scattering is far more complicated than electron 
scattering ||, as excitations and ionization of Ps are expected to play a significant role for each 
target states. Besides if the target is vulnerable to excitations, prediction of total cross section 
becomes extremely difficult as the number of scattering channels grow as N 2 , where N excited 
states of both Ps and target are considered. 

Due to internal charge and mass symmetry, Ps atom yields zero elastic and even-parity tran- 
sition potentials in the direct channel and Ps scattering is dominated mainly by exchange corre- 
lation ]7|, |8], [|. This eventually complicates the convergence of a conventional coupled- channel 
description [[/] . For Ps scattering with molecular target we have additional complication due to 
the presence of the nuclear degrees of freedom. The presence of three charge centers complicates 
the solution scheme of the coupled-channel equations. The degree of complication can be real- 
ized from the mathematical analysis of ref. [|K| for Ps formation in positron- hydrogen-molecule 
scattering. Because of mathematical complication, this approach has not been further pursued 
in numerical analysis and there exists no successful calculation of Ps formation in positron- 
hydrogen-molecule scattering. We recall that the mathematical complication of a dynamical 
coupled-channel study of Ps formation in positron-hydrogen-molecule scattering is similar to 
that of Ps-H 2 scattering. 

Here we undertake a theoretical study of Ps scattering by H 2 . We employ a hybrid approach 
of treating the molecular orientational dependence as parameters in the coupled equations and 
then perform the partial- wave expansion fTI| . The resulting one- dimensional equations are 
solved at the equilibrium nuclear separation by the matrix-inversion technique and the partial- 
wave cross sections are numerically averaged over molecular orientations. The present approach 
makes the Ps-H 2 scattering problem easily tractable. 

Recently, we suggested a regularized electron-exchange model potential |], || and demon- 
strated its effectiveness in exchange-dominated Ps scattering by performing quantum coupled- 
channel calculations using the ab initio framework of close- coupling method for both simple 
and complex targets |7], 12, 13, 14, 15]. In our initial calculations we used a nonsymmetric 
model-exchange potential for Ps scattering by H |T2||, He M and H 2 R O] and obtained rea- 



sonably good agreement with experiment on He and H 2 . In a subsequent application of Ps 
scattering by H, it was found that a time-reversal symmetric form of the exchange potential 
leads to far superior result than the nonsymmetric form both in qualitative and quantitative 
agreement with accurate variational calculations on H 0. The symmetric potential also led 
to very good results for low-energy cross sections for Ps scattering by He, Ne, and Ar in 
excellent agreement with experiment Q. 

In view of the above we reinvestigate the problem of Ps scattering by H 2 using the symmetric 
exchange potential employing the three-Ps-state [Ps(ls,2s,2p)] coupled-channel model for elastic 
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and Ps(2s,2p) excitations. We solve the coupled-channel equations by the above scheme and 
report partial cross sections for Ps(ls,2s,2p) excitations. We also calculate cross sections for 
Ps excitations to 6 > n > 3 states and Ps ionization using the first Born approximation. A 
target-elastic total cross section is calculated by adding the above partial cross sections. We use 
the present symmetric exchange potential to calculate the target-inelastic Born cross sections 
for excitation to B X E+ and b 3 E+ states of H 2 by Ps impact. We present a total Ps-H 2 cross 
section by adding the above target-elastic and target-inelastic results. 

In section 2 we present the theoretical formulation in the body-frame fixed-nuclei approx- 
imation. In section 3 we present the numerical results. Finally, in section 4 we present a 
summary of our findings. 



2 Theoretical Formulation 

The total wave function \1/ of the Ps-H 2 system is expanded in terms of the Ps and H 2 quantum 

states as 

*(r ,r 1 ,r 2 ,x;R) = ^ab(so)y?a(t )?Mri, r 2 ; R 

a, 6 

where s = (x + r )/2, t = (x — r ), x (r ) is the coordinate of the positron (electron) of the 
Ps atom, ri and r 2 are the coordinates of the electrons of H 2 , 2R is the internuclear separation 
of H 2 , tp a the wave function of Ps, ipf, the wave function of H 2 , with b denoting the electronic 
configuration of H 2 and a denoting the quantum state of Ps. Here A denotes antisymmetrization 
with respect to Ps- and target-electron coordinates and F a t, is the continuum orbital of Ps with 
respect to the target. The spin of the positron is conserved in this process and the exchange 



profile of the Ps-target system is analogous to the corresponding electron-target system |L4 . 

Primary complication arises in the coupled-channel study in retaining both the summations 
over Ps and target states in the coupling scheme. As we shall see later that the contribution 
of the individual target-inelastic channels to cross section is one order small compared to that 
of target-elastic channels. So we exclude the summation over target states from the coupling 
scheme and treat cross section for target excitations separately using Born approximation. In 
this work we are mostly interested to see whether the model can accommodate the measured 
total cross sections of Ps-H 2 scattering. So we retain the Ps(ls,2s,2p) excitations in the coupled- 
channel model and treat the remaining excitations by Born approximation. The first Born 
calculation with present regularized exchange leads to results reasonably close to coupled- 
channel calculation at medium energies. Hence it is expected that this calculation should yield 
a fairly good picture of the total cross section except near the Ps-excitation thresholds. 

Projecting the resultant Schrodinger equation on the final Ps and target states and averaging 
over spin, the resulting momentum-space Lippmann-Schwinger scattering equation in the body- 
frame representation for a particular total electronic spin state S can be written as 

fijy, k; R) = , k; R, - JL £ / ^ S^Ml^m (2) 

where f^, a are the scattering amplitudes, and £>f/ a the corresponding Born potentials. k a " is the 
on-shell momentum of Ps in the intermediate channel a". We use units h = m = 1 where m 
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is the electronic mass. For PS-H2 target-elastic scattering there is only one scattering equation 
(TJ) corresponding to total electronic spin S = 1/2. The input potential of @ is given by 

B'Jtik', k; R) = B°jy, k; R) - B* a (k', k; R) (3) 

where B D is the direct Born potential and B E is the exchange potential. As in the ground 
electronic state X l H+ of H 2 , the total electronic spin is zero, there will contribution of one 
target electron to B% a (]t', k; R) Q. 

For the electronic ground state X 1 S+ of H2 we use the wave function of the form ^,(ri, r 2 ; R) 
= laf\l)laf\2) = A^ 2 f/ (r 1 ;R)f/ (r 2 ;R), where N = [2(1 + T)]~ 1/2 with T the overlap 
integral ||| and C/ (r; R) = (5 3 /vr) 1 / 2 [ eX p(-5|r - R|)+exp(-<5|r + R|)] with 5 = 1.166 ||, ||. 
For Ps we use the exact wave functions, e.g., for Is state cp(r) = exp(— /3r)/v / 87r with (3 = 0.5. 
For the excited states of H 2 we take the configurations i/v(ri, r 2 ; R) = [1<t£ 6 '(1)1<t4 6 '(2) ± 
lcr^(2)lcr^ 6 ^(1)], where + (— ) corresponds to the spin singlet (triplet) state. All wave functions 
of H 2 are taken from ref. [ |T9|j . 

The direct Born potential for Ps transition from state a to a' and H 2 transition from b to b' 



can be rewritten in the following convenient factorized form [^0 



£ a V a6 (k f ,k i; R) = ^ f dt&(t)[e^ - e^Wfrait) 



Q 2 

J dr 1 dr 2 ^(ri,r 2 ;R) 



2cos(Q.R) - ^ e iQ - r " 



71=1 



V 6 (r l5 r 2 ;R). (4) 



Next we describe the electron-exchange model potential for Ps-H 2 scattering. We develop 
the model exchange potential from the following term P, |T^] : 

£^ ab (k f ,k i; R) = — [ dxdT Q d ri dT 2 e-^ x+r ^V a ^ ~ ri)#(r , r 2 ; R) 1 

7r J \r — ri| 

x ^(ri,r 2 ;R)v9 a (x-ro)e tki -( x+r °)/ 2 . (5) 

After removing the nonorthogonality of the initial and final wave functions of (|5|) and some 
straightforward simplification the model exchange potential for a general target-elastic transi- 
tion becomes || |§ 

e„V.( k( , k ,;R) = ^ - m + 4( -g"; g - K)j2 J ^W^V.W 

I dr,dr i,;(r , r 2 ; R)e ,Q "V t (r„, r 2 ; R) (6) 

(^ + ^)/8 + 4 ^r +/g+ ^)/2 /^^ ,/a ^ 

I dr 2 laf (r 2 )laf (r 2 ) / dr laf (r K Q - r °laf(r ) (7) 



x 



x 



where Z and /' are the angular momenta of the initial and final Ps states and Q = ki — kf . 
Here 8^ and j3 a are the parameters of the H 2 and Ps wave functions in the initial state. The 
parameter /3 2 / corresponds to the final-state binding energy of the Ps atom and is taken as 
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zero while considering exchange for the Ps ionization channel. For target-inelastic processes, 
following the prescription outlined in refs. || U the exchange potential takes the form 



In (|3|) the additional indices 6 and 6' are introduced on the molecular orbitals la g and la u 
to distinguish the initial and final states. The model exchange potentials (|7]) and (@) may 
be considered as a generalization of a similar potential suggested by Rudge pi] for electron- 



atom scattering. The nonlocal and time-reversal symmetric exchange potential (|7|) has a very 
convenient form and is expressed as a product of form factors of Ps and H 2 . Both the direct and 
exchange amplitudes have been factored out in terms of Ps and target "form-factors" leading 
to a substantial simplification of the theoretical calculation. In our previous study of Ps-H 2 
scattering || |H| the prefactor of the exchange potential was not time-reversal symmetric. In Eq. 
(0) we have restored time-reversal symmetry as in ref. || which is found to provide significant 
improvement in the results. Although, for electron-molecule scattering several model potentials 



are found in the literature |22|, p3|l , there is no other convenient model exchange potential for 
Ps scattering. 



3 Numerical Procedure and Results 

In the body-frame calculation the coupled-channel equations are solved at the equilibrium 
nuclear separation 2Rq = IAclq. The polar and azimuthal angles Or and <J)r of R are taken as 
parameters in the coupled equations O]. This reduces (|2|) to the following form 

fa'a (k,k)-B Q , a (k,k)-— ^Jdk k 2 /i/A _ k ,2 /A + i0 • I 9 ) 

After standard partial-wave projection the three-dimensional coupled scattering equations (|9|) 
are first reduced to coupled one- dimensional integral equations in momentum space. The one- 
dimensional equations are then discretized by Gauss-Legendre quadrature rule and solved by 
the matrix inversion technique. A maximum of forty points are used in the discretization of 
the integrals. The discretized coupled-channel equations are solved for eight to ten discrete 
values each of polar and azimuthal angles Or and <pR. This leads to a convergence of cross 
sections up to three significant figures. For targets with more charge asymmetry and for polar 
molecules we expect that more points will be required for angular averaging. The present 
averaging amounts to solving the coupled set of scattering equations sixty four to hundred times. 
Although this procedure increases the computational (CPU) time, mathematical complications 
of tedious (and untractable) angular-momentum analysis |TtJ are thus replaced by a tractable 
and convenient calculational scheme. Finally, the partial-wave cross sections are numerically 
averaged over molecular orientation using Gauss-Legendre quadrature points for both polar and 
azimuthal angles Or and <pR. Maximum number of partial waves included in the calculation 
is 12. Contribution of higher partial waves to cross section is included by corresponding Born 
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terms. These Born cross sections are also numerically averaged over all molecular orientations 
in a similar fashion. 

In figure 1, we show the angle-integrated target-elastic cross sections for elastic, Ps(2s+2p), 
Ps(3 < n < 6) excitations and ionization of Ps. As expected and observed in previous calcula- 
tions |7|, ||, [12], [L3| with other targets, the contribution of the Ps ionization channel to the cross 
section plays a dominant role from medium to high energies as can be seen in figure 1. The 
detailed angle-integrated partial cross sections of the Born and three-Ps-state calculation are 
tabulated in table 1. Near 10 eV, 20 eV and 30 eV the total Born cross sections for Ps(ls,2s,2p) 
are found to be nonconvergent by 28%, 9% and 4%, respectively. Near 10 eV, 20 eV and 30 eV 
the total Born cross sections for Ps(n = 2) excitations are found to be nonconvergent by 18%, 
6% and 3%, respectively. At 60 eV the Born and three-Ps-state results are essentially identical. 
The nonconvergence of Ps(n > 3) excitations calculated using Born approximation are expected 
to lie within the limit set by Ps(n = 2) cross sections above. It is expected from experience 
that the ionization cross section calculated using Coulomb Born will be more converged than 
the Ps(n = 2) Born cross sections. 

Table 1: PS-H2 partial cross sections in units of ita^ at different positronium energies using 
the Born approximation and three-Ps-state calculation 



E (eV) 


Ps(ls) 


Ps(ls) 


Ps(2s) 


Ps(2s) 


Ps(2p) 


Ps(2p) 


Ps(n > 3) 


Ps-ion 




Born 


3-St 


Born 


3-St 


Born 


3-St 


Born 


Born 


0.068 


23.72 


3.79 














0.612 


17.56 


3.16 














1.45 


11.84 


2.47 














3 


6.71 


1.60 














4 


5.03 


1.17 














5 


3.94 


0.81 














6 


3.18 


1.01 


0.098 


0.23 


3.46 


2.28 








7 


2.63 


1.05 


0.101 


0.18 


3.60 


2.52 


1.28 


0.11 


8 


2.22 


1.03 


0.092 


0.14 


3.36 


2.51 


1.43 


1.15 


10 


1.65 


0.94 


0.072 


0.090 


2.80 


2.27 


1.29 


3.03 


12.5 


1.21 


0.81 


0.052 


0.068 


2.24 


1.93 


1.05 


4.34 


15 


0.92 


0.68 


0.041 


0.054 


1.85 


1.65 


0.87 


4.94 


20 


0.59 


0.49 


0.026 


0.035 


1.34 


1.25 


0.63 


5.18 


25 


0.40 


0.36 


0.017 


0.023 


1.04 


0.99 


0.49 


4.93 


30 


0.29 


0.27 


0.012 


0.017 


0.84 


0.81 


0.40 


4.55 


40 


0.17 


0.16 


0.007 


0.009 


0.60 


0.58 


0.28 


3.80 


60 


0.07 


0.07 


0.003 


0.004 


0.37 


0.37 


0.17 


2.72 



Now we concentrate on some target inelastic processes. Each target inelastic transition is 
accompanied by elastic, excitation and ionization of Ps and hence by an infinite number of 
possibilities. Here to account for target inelastic processes we consider Ps(n = 1 — > 6) discrete 
excitations and ionization of Ps using the first Born approximation. Contribution of higher 
discrete excitations of Ps are expected to be insignificant and are neglected in the calculation. 
We calculate the cross sections for the inelastic transition Ps(ls) + H 2 (X — > Ps* + H 2 (B 
X S+) where Ps* represents the ground and the excited states of the Ps atom. These Born cross 
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sections are also calculated at a equilibrium internuclear separation 2Rq = 1.4ao and finally 
averaged over angular orientations of the target. In figure 2 we display the Born contribution 
of partial cross sections for transition to different Ps states while the hydrogen molecule is 
excited to the H 2 (B state. The total cross section summing the different contributions 

is also shown for this target-inelastic process. When compared with the corresponding total 
cross section calculated with the nonsymmetric exchange potential (not shown in figure 2) 
we find marginal change at low energies (below 20 eV) and basically no change at medium to 
high energies. This is quite expected as the effect of exchange dies out at higher energies and 
these cross sections are controlled by the direct Born potentials. 

Next we calculate the cross sections for the inelastic transition Ps(ls) + H 2 (X — > Ps* 
+ H 2 (b 3 £„) using the first Born approximation, where Ps* represents the ground and the 
abovementioned excited states of the Ps atom. These cross sections are again calculated for 
equilibrium internuclear separation 2R = lAa Q and averaged over molecular orientations. In 
figure 3 we display the contribution of partial cross sections for transition to different Ps states 
while the hydrogen molecule is excited to the H 2 (b 3 S^) state. The total cross section summing 
the different contributions is also shown for this target-inelastic process. There is significant 
change when we compare this total cross section calculated with the symmetric exchange po- 
tential with our previous result [0] obtained with the nonsymmetric exchange potential. This is 
quite expected as this process is purely exchange dominated. The qualitative and quantitative 
differences between the two total cross sections reveal the importance of using the symmetric 
exchange potential. 

In figure 4, we exhibit the total cross section obtained from target-elastic processes (Ps 
excitation up to n = 6 and Ps ionization) plus target excitations to B : £+ and b 3 £+ states 
H [13| and compare with the total cross section measurements of ortho-Ps scattering from H 2 
|2], |[]. The agreement with the experimental cross sections is highly encouraging. Here, we have 
considered only two lowest target-inelastic processes and their combined effect on the total cross 
section. This gives an indication that the inclusion of remaining important target-inelastic cross 
sections might give a better agreement with measurement. 

However, the present theoretical peak in total cross section is shifted to a lower energy 
compared to experiment. This trend was also found in our calculation of Ps-He scattering 0. 
In the energy range of 7 to 15 eV, the total cross section overestimates the measured data. 
This is due to the fact that Ps(n > 3) excitations and ionization have been treated in the 
first Born approximation framework. The neglect of other target-inelastic channels and the 
first Born calculation for higher excitations and, in particular, ionization of Ps are supposed 
to be responsible for the shift in the theoretical peak. A dynamical calculation for higher Ps 
excitations and ionization should reduce the theoretical total cross section in the intermediate 
energy region, while the inclusion of further target-inelastic channels will increase the cross 
section for medium to higher energies. These two effects are expected to shift the peak in the 
total cross section to higher energies and lead to a better agreement with experiment. 

A previous first Born description of Ps-H 2 scattering [24| with Born-Oppenheimer exchange 



25] led to unphysically large cross section at low energies. The poor performance of that 
scheme in Ps scattering compared to electron scattering is due to the fact that in the absence 
of the direct potential (zero for elastic scattering), the Born-Oppenheimer exchange potential 
solely determines the cross section. This clearly shows the very unrealistic nature of the Born- 
Oppenheimer exchange potential at low energies. 
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In previous studies |7], |8|, [L4| of Ps scattering we introduced a parameter C in the exchange 



potential for obtaining a more quantitative fit with experiment essentially by replacing the 
denominator term in (0) by (kj + kf)/8 + C{8 2 + 8%, + (3% + /3y)/2. In the original form the 
constant C = 1 and we have used this value in this study. However, it can be varied slightly from 
unity to obtain a precise fit of a low-energy scattering observable (experimental or variational), 
as has been done in some applications of model potentials [p3|| . A variation of C from unity 
leads to a variation of the average values for square of momenta |8], |T2]| , which were taken 
as the binding energy parameters (5 2 ,/3 2 etc.) in the expression for the denominator in (0). 
This variation, in turn, tunes the strength of the exchange potential at low energies. At high 
energies this model potential is insensitive to this parametrization and leads to the well-known 
Born-Oppenheimer form of exchange [p5] . 



4 Summary 

To summarize, we have used a time-reversal symmetric form of the nonlocal model potential for 
exchange and applied it to the study of Ps-H 2 scattering. We have also presented a simplified 
prescription for performing coupled- channel dynamical scattering calculation with molecular 
target using a body-frame fixed-nuclei scheme. With this prescription we have performed a 
three-Ps-state coupled-channel calculation of target-elastic Ps-H 2 scattering. Higher excitations 
and ionization of the Ps atom are treated using the first Born approximation with a regularized 
exchange. We also calculated cross sections for two target-inelastic excitations of H 2 (B X E+ 
and b 3 S+) using first Born approximation with present exchange considering Ps excitations 
(n = 1,...,6) and ionization. Considering the fact that we have considered only two target- 
inelastic processes, the present total cross section is in encouraging agreement with experiment. 
The tractability of of the present dynamical calculational scheme for molecular targets and the 
success of the present time-reversal symmetric electron-exchange potential in describing Ps-H 2 
scattering should stimulate further investigation with both. 

We thank the Conselho Nacional de Desenvolvimento Cienti'fico e Tecnologico, Fundagao de 
Amparo a Pesquisa do Estado de Sao Paulo, and Financiadora de Estudos e Projetos of Brazil 
for partial financial support. 
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Figure Caption: 

1. Angle integrated target-elastic PS-H2 partial cross sections at different positronium ener- 
gies: elastic (solid line) and Ps(2s+2p) excitation (dashed-dotted line) from three-state coupled- 
channel model at low energies interpolated to exchange Born at high energies, Ps(6 > n > 3) 
excitation (dashed-double-dotted line), Ps ionization (dashed line) from exchange Born. 

2. Angle integrated target-inelastic PS-H2 partial and total cross sections at different positro- 
nium energies to H 2 (B *£+) state using the first Born approximation: Ps(ls) (dashed-dotted 
line), Ps(2s+2p) (dashed-double-dotted line), Ps(6 > n > 2) excitation (dotted line), Ps ion- 
ization (dashed line), and total (full line). 

3. Angle integrated target-inelastic PS-H2 partial and total cross sections at different positro- 
nium energies to H 2 (b 3 S+) state using the first Born approximation: Ps(ls) (dashed-dotted 
line), Ps(2s+2p) (dashed-double-dotted line), Ps(6 > n > 2) excitation (dotted line), Ps ion- 
ization (dashed line), total (full line), total cross section with the nonsymmetric potential from 
ref. || (dashed-triple-dotted line). 

4. Total Ps-H 2 cross section at different positronium energies: total target elastic cross 
section from figure 1 (dashed line), total target-elastic (from figure 1) plus target-inelastic 
(from figures 2 and 3) cross section (solid line), experimental data (solid circles from ref. ||, 
solid square from ref. [|J). 
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Figure 1 
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